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Abstract 

We resum to next-to-leading order the distribution in the ratio of the invariant hadron mass mx to the total 
hadron energy Ex and the distribution in mx in the semileptonic decays B — > X u lv. By expanding our 
formulas, we obtain the coefficients of all the infrared logarithms at 0{a 2 s ) and of the leading ones at 0(a|). 
We explicitly show that the relation between these semileptonic spectra and the photon spectrum in the radiative 
decay B — » X s ^ is not a purely short-distance one. There are long-distance effects in the semileptonic spectra 
which are not completely factorized by the structure function as measured in the radiative decay and have to 
be modelled in some way. 
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1 Introduction 



Semileptonic B decays 

B -» X u + I + v, (1) 

where X u is any hadron final state coming from the fragmentation of the up quark, are interesting processes 
for the study of strong interactions as well as of weak interactions. The computation of the spectra in Q is 
often non trivial because of the presence of double infrared logarithms in the perturbative expansion, which 
formally diverge in the endpoints and therefore must be resummed to all orders in the QCD coupling as pQ. 
In all generality, large logarithms come from the so-called threshold region, defined as the one having 

ra x « E x < m b , (2) 

where nix and Ex are the invariant mass and total energy of the final hadron state X u and mb is the beauty 
mass. We find useful to summarize here the main results of [21 EJ- The infrared logarithms in process can 
be organized in a series of the form 

oc 2n 

(3) 



E[«; a(Q)\ = 1 + ££s nfc a n (Q) logM -) 

n=lfc=l ^ U ' 



and can be factorized into the universal QCD form factor X. The E^'s are numerical coefficients whose explicit 
expressions can be obtained from and a = as is the QCD coupling. Q is the hard scale of the process and 
is determined by the final hadron energy: 

Q = 2E X . (4) 

We have defined the hadron variable 



l-Vl-(2r»^/Q) 2 fm x \ 2 ... <n (K . 

u= -, = ~ — (0<u<f), (5) 

1 + yj 1 - (2m x /Q) 2 V Q J 

involving the ratio of the invariant hadron mass to the hard scale, where in the last member we have taken the 
leading term in the threshold region (|2J) only. As it is well known, the QCD form factor in eq. © contains at 
most two logarithms for each power of a, coming from the overlap of the soft and the collinear region in each 
emission. Note that the hard scale Q enters the argument of the infrared logarithms 1/u ~ Q 2 /m x as well as 
the argument of the running coupling a = a(Q). One can obtain a factorized form for the triple differential 
distribution — which is the most general distribution in process Q — from which all other spectra can be 
obtained by integration: 

1 f u d 3 T 

du = C [x, w; a(w nib)} S [u; a(w mb)] + D [x, u, w\ a(w nib)] , (6) 



T J dxdwdu 

where 



w = — (0<w<2); x = — (0<x<l). (7) 

nib m b 

r is the total semileptonic width, C [x, w; a] is a short-distance coefficient function independent on u and 
D[x,u,w;a] is a remainder function not containing infrared logarithms (i.e. short-distance dominated) and 
vanishing for u — > as well as for a — > 0. The explicit expressions of these functions have been given in py. 

The properties of semileptonic decay spectra are best understood comparing them with the simpler radiative 
decay 

B - X s + 7 . (8) 

In such decay we have indeed: 

Q = m b (l-JL + = m b ( 1 + ~ m b , (9) 

\ ml ml) \ ml) 
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where is the 4-momentum of the real photon — in general of the probe. In the radiative decays (JSJ the hard 
scale is therefore independent on the kinematics and is fixed by the beauty mass. In the semileptonic decay Q, 
q^ is the dilepton momentum and we have the more general situation < q 2 < m 2 ; the hard scale is given in 
this case by 

Q~m h (\-?—\ (10) 



\ m 6/ 

and depends on the dilepton invariant mass squared q 2 : it cannot be identified with the heavy flavor mass m b . 
Kinematic configurations with 

m x < Q ~ m b (11) 

as well as with 

m x < Q < m b (12) 
are possible. In fact, in the radiative decays, for example, the average hadron energy is 

(Ex)rd = \m b [1 + 0(a)], (13) 

while in the semileptonic ones we have the smaller value [3] 

(Ex)si = ^m h [1 + 0(a)}. (14) 

In general, semileptonic spectra may or may not involve integration over the hadron energy Ex and, according 
to 001], have a different infrared structure in the two cases. In we have studied in detail the simpler case of 
the distributions not integrated over the hadron energy, i.e. over the hard scale Q, which have the same infrared 
structure of the hadron invariant mass distribution in the radiative decay (JHJ: 



1 r ts rIV / oo 2n \ 

f" Idf dt ' s = Cli{a) 1 + E E "( m ")" l0 S fe T + ^«(*-5 «)• ( 15 ) 

R JO s \ n=lk=l S / 

We have defined t s = m 2 x /m 2 = 1 — x 7 , with .t 7 = 2E 1 /m b , Tr is the total radiative width, Cji(a) is a 
short-distance coefficient function and Dn(t s ;a) is a short-distance remainder function. 

In this paper, we attack the distributions integrated over the hadron energy, which have a more complicated 
infrared structure than the one in JSJ|. 

In sec. we resum to next-to- leading order (NLO) the distribution in the hadron variable u defined before. 
The infrared logarithms appearing in the perturbative expansion of this spectrum, 



F / tT 7 ^' = Cu ^ [ 1 + ^2^2^Unk a (m b ) n \og k - J + Du(u;a), 

J® U \ n=lfc=l u ) 



(16) 



coincide in first order with those in the decay JSJ: £[/i2 = S12 and Sjjh = Sn, while they differ in higher 
orders. Cu{a) is a short-distance coefficient function and Djj(u; a) is a short-distance remainder function whose 
explicit expressions will be given in sec. (01- 

In sec. we compute to NLO the distribution in the variable 

TTl 2 

t = Hr (o<t<i), (17) 

i.e. the distribution in the invariant hadron mass squared. The infrared logarithms appearing in this distribution, 

lk a(m b ) n log fc i ] + D T (t;a), (18) 



J0 \ n=lk=l 



2 



differ also at the 0(a) single logarithm level from the corresponding ones in (JSJ: Stii 7^ S n . We also define 
a non-minimal factorization-resummation scheme which seems to have better convergence properties of the 
perturbative scries than the minimal one. 

We define for both spectra effective form factors which resum the large logarithmic corrections. These form 
factors involve a convolution of a process-dependent coefficient function with the universal QCD form factor £ 
entering the radiative decay J5J) and the triple differential distribution in the decay . The definition of the 
effective form factors we give is, in a sense, non perturbative, because we do not look at the explicit perturbative 
expansion of X. 

Finally, in sec. ^ we draw our conclusions and we consider generalizations of our results. 



2 Distribution in the hadron mass/energy ratio 

In this section we compute the resummed distribution in the variable u (defined in the introduction) to next- 
to-leading order (NLO). That is accomplished by integrating the resummed double distribution in u and w 
obtained in sec. (4) of |3j over w: 

ldT _ f 1+u 1 dPT 
r du J Q r dwdu 



By replacing the resummed expression on the r.h.s. of eq. I|19|l . we obtain: 



1 dT 

r du 



1+u 



1+u 



dw Ch(w; a) a [u; a(w mt,)] + / dw dniu.w] a) 
'o Jo 

Ch(w; a) is a short-distance coefficient function having an expansion in powers of a: 

C H {w; a) = C { °\w) + aC [ ^{w) + a 2 [w) + 0{a 3 ), 

with 



(20) 



(21) 



4 0) H - 



2w 2 (3-2w); 
— 2w 2 (3-2w) 



-r • / \ , , M \ 35 9-410, 
Li 2 {w) + logiolog(l - 10 J 6 - 4w W 



(22) 
(23) 



where Cf = {N% — 1)/(2N C ), N c = 3 is the number of colors and a = a(mt,). a (u; a) is the differential QCD 
form factor: 

a(u; a) — — a), (24) 
du 

with a) being the cumulative form factor considered in the introduction, 

E(w; a) = 1 + qS (1 >(«) + a 2 E (2) (w) + 0(a 3 ) 

and 



7T 



-log u + -logtt 



(25) 



(26) 



Finally, djj(tt, w; a) is a short-distance remainder function whose explicit expression is not needed here. We are 
interested in the threshold region which can also be defined as the one having 



u < 1. 



(27) 



Since large logarithms originate only from the first integral on the r.h.s. of eq. (|2U[l. in order to isolate them, 
let us neglect at first the contribution from the remainder function du(u,w; a). The small terms for u — * 
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will be included later on, by expanding the resummed expression and comparing with fixed-order spectrum, 
as discussed at the end of sec. 2 of 3 . The integration of the first term on the r.h.s. of eq. ((2U(I over the 
hadron energy is not trivial because the coefficient function Ch(w; a) depends on w as well as the QCD form 
factor a[u; a(umj)], which depends on w via the scale of the coupling a — a{wm,b). Unlike the distributions 
considered in we cannot factor out in this case the universal form factor a. 

Since it is technically simpler to deal with partially-integrated form factors rather than with differential 
ones, let us define the event fraction 

f u 1 dT 

Ru(u) = j o f^jdu', (28) 

which has the end-point values: 

Ru(0) = 0; Ru(l) = 1. (29) 

The spectrum in u is trivially obtained by differentiation: 

1 dT d 

fTu-Tu R ^- W 

Integrating both sides of eq. 1)20(1 . we obtain: 

Ru [u; a{m%)\ — I du 1 I dwCH(w;a)a[u';a(wmb)] + 0(u,a) 
Jo Jo 

du' i dwCH(w;a)a[u';a(wm,b)]+ / du' / dw Ch{w; a) a[u'; a(w mb)] + 
o Jo Jo Ji 

+ 0(u;a), (31) 

where by 0(u; a) we denote terms which vanish for u — > as well as for a — > 0. The second integral in the last 
member of the r.h.s. extends to a kinematic region in w which is 0(u) and therefore can be dropped in the 
threshold region; we can therefore assume tree-level kinematics: 

< w < 1. (32) 

By exchanging the order of the integrations in the remaining integral, we obtain: 

Ru[u;a(mb)] = / dw Ch(w; a) E[u; a(w mb)] + 0(u;a). (33) 
Jo 

Substituting the explicit expressions for the coefficient function given in eq. (|2f I) and of the QCD form factor 
given in eq. ((25(1 . we obtain: 



dwCH(w;a)T,(u; a) = 1 — — 



1 2 7, 335 
-log u + -logu + — 



+ 0{a 2 ). (34) 



The next step is to factorize the event fraction Rjj{u; a) into: 



• a QCD form factor S[/(u; a) containing the long-distance contributions, i.e. the logl/w terms diverging 
for u — > 0; 

• a coefficient function Cu{ct) containing the constant terms for u — > 0; 

• a remainder function Du(u; a), collecting the left-over small contributions O{u;o), vanishing for u — > 
and for a — > 0. 
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Let us write therefore: 

Ru (u; a) = Cu{ot) S[/(u; a) + Djj(u, a). (35) 
The coefficient function, the form factor and the remainder function can all be expanded in powers of a: 



Cu{a) = 1 + adfi + a 2 C ( jj ] + 0(a 3 ); (36) 
Eu(u; a) = 1 + aE^u) + a 2 ^\u) + 0{a 3 ); (37) 
Du(u; a) = aD$\u) + a 2 D (2) {u) + 0(a 3 ). (38) 



The knowledge of soft-gluon dynamics allows the resummation of the dominant terms to all orders in a in Ej/; 
this is instead not possible for the coefficient function and the remainder function, which are not long-distance 
dominated and for them one has to use truncated expansions. 

The above conditions do not completely specify the form of the form factor, of the coefficient function and of 
the remainder function, so we have to select a factorization scheme. Let us choose a minimal scheme, in which 
the form factor contains only logarithmic terms: 



Et,(u; a) = 1 + J2 ^ Unk ° n LK ( 39 ) 

where 



n=l fc=l 



L = log-. (40) 
u 

From the above definition, it follows that has the same normalization as E: 

a) = 1, (41) 
which holds to any order in a. We obtain in first order: 



7T 



-log U + -l0gM 



(42) 



(i) 335 C F 

° u 144 V (43) 



Note that 



££>(«) = E«(«). (44) 



As we are going to explicitly show later in this section, this property does not hold in higher orders of a. For 
a(m&) = 0.22, we have a first order correction to the coefficient function of — 21.7%. 

The remainder function is obtained by imposing consistency between the resummed expression and the 
fixed-order one (matching). By expanding to first order in a the resummed expression (|35|l and imposing the 
equality with the full 0(a) result, one obtains 0]: 



C 

D u\ u ) = ~~ — ^TTTTT— 1 ' — 1 ' ~ — — oi7 — ul °S u ■ ( 45 ) 



7T 



u(15624 - 2688u - 1352u 2 + 141u 3 ) 21 - 84u - 29u 2 + 6u 



3 



5040 210 



As required, the remainder function goes to zero in the elastic point u — (as ulogu). Taking u = 1 in 
eq. Ij35(l . we obtain the following relation between the coefficient function and the remainder function in the 
upper endpoint: 

Cu(a) = 1 - D v (l;a), (46) 
which holds to any order in a and is verified in first order (see eqs. I|43|l and (14511 1. 
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The minimal scheme has been defined above by looking at the explicit form of the event fraction Rjj(u] a) 
as a power series in a: one reorganizes the series picking up the logarithmic terms and putting them into the 
effective form factor £(/. It is also possible to give a different definition of the minimal scheme which does not 
make use of the explicit expansion of the event fraction. Since 

l 

dw Ch{w; a) £ [u; a(wrrib)} = Cu(pt) Y,u{u\ a) (47) 

o 

and 

Sc/(1; a) = 1, (48) 
we can define the coefficient function to all orders as: 

Cu(a) = [ dwC H (w; a). (49) 



o 



The effective form factor is therefore written as: 



J dwC H (w; a)T,[u; a(wm b )} 

Eu (u; a) = -a . (50) 

J dwC H {w] a) 

The definition (|50|l of the minimal scheme has the following phenomenological advantage. One may wish to use 
for £(u;a), instead of the perturbative expression, for example the result of a fit to some experimental data 
or a non-perturbative model 4 . In these cases, Y,{u\w) does not depend on the coupling and therefore cannot 
be expanded in a, but the effective form factor £[/ in the minimal scheme can still be computed by means of 
eq. (53 . 



The representation of the effective form factor £j/ given in eq. (|50|l allows us to make a few general comments: 

• £[/(m; a) factorizes all the threshold logarithms in the spectrum but is, unlike £(u), process dependent. 
That is because it involves the convolution of the universal form factor E(u; a) with the process-dependent 
coefficient function over all the hadron energies. Ch(w; a) has the role of a probability distribution: it 
gives the probability for the hadronic subprocess with hard scale Q — w m b to occur. In the u distribution 
hadronic subprocesses with all the possible hard scales Q from zero up to m b do contribute, while in the 
radiative decay © Q is kinematically fixed to the upper value nib- The relation between these two spectra 
therefore is not a purely short-distance one: to relate these two distributions one has to model in some 
way the variation of the form factor cx(Q)\ with Q ranging from nib down to zero. In agreement with 
physical intuition, the problem in the computation of the u distribution is the estimate of the contributions 
from small hard scales, Q <C nib, where perturbation theory is expected to fail. However, since Ch oc Q 2 , 
small hard scales give a small contribution to the total. We may say that the u spectrum is "protected" 
from long-distance effects related to large logarithms with a large coupling, i.e. related to the region in 
eq. (O; 

• since a(wrrib) — a(rrib) + 0(a 2 ), the effective form factor and the universal form factor coincide in first 
order, as already found by explicit computation (see eq. I|44ll ): 

• were not for the dependence of the form factor £ = £[ii; a(wrrib)] on w through the running coupling, 
a = a(wnib), the effective form factor Ej/ would coincide with the universal one S to all orders. 

Let us now discuss the higher orders in the perturbative expansion of Tiu(u; a). One has to insert in eq. (|50|l : 

• the QCD form factor S, whose next-to-next-to-leading order corrections (NNLO) are now well established 
due to the recent revaluation of the resummation constant Di in 0; 



4 For the coefficient function, one can still use the perturbative result. That is because the coefficient function, unlike the form 
factor, is short-distance dominated, and therefore its perturbative evaluation is more reliable. 



(i 



• the coefficient function Cjj(w; a), which is known to 0(a), i.e. to NLO (cfr. eq. iJHJ), 



and perform the integration over w. To obtain the truncated expansion of £[/, one simply replaces the truncated 
expansions of £ and of Cu , expands the product and integrates term by term. 

The exponential structure of the threshold logarithms in £(u; a) is partially spoiled in Y,jj(u; a) because 
of the integration over the hard scale Q = wrrib, but it is not completely ruined. In order to simplify the 
representation of the large logarithms, it is therefore convenient to introduce the exponent of the form factor 
also in the effective case, i.e. to define Gjj as: 

Ep = e Gu . (51) 

Expanding the exponent in powers of a up to third order included, we have an expansion of the same form of 
G, defined in [3], that is: 

oo n+1 

G v {u; a) = ^ Y GlJnk a " L * = G ui2aL 2 + G un aL + G U23 a 2 L 3 + G U22 a 2 L 2 + G U2 ia 2 L + ■■■ , (52) 

n=l k=l 

where 

L = log- > (53) 
u 

and a = a(m&). In order to explicitly see the effects of the integration over the hard scale, let us give the 
coefficients Guij up to the third order in terms of the corresponding coefficients Gij of the radiative decay JHJ 
and of the resummation constants Ai , Bi , Di and of the /3-function coefficients f3i |2j : 





= G 




= G 


G[/23 


= G- 


Gu22 


= G; 


G\J2\ 


= G; 


Gu34 


= G : 


Gu33 


= G; 


G(J32 


= G; 



(54) 
(55) 
(56) 

y^i/3o; (57) 
^o(Si + £>i); (58) 

5 
6 

5 1 A ( nn ,t 9 _ _ ,t \ 5 ,^ 9 / ^ ,~» ^ \ 5 
-A 2 3 

6 H 36 



(59) 

Aif3 2 ; (60) 
^ A t (23 13 2 + 15 0t) - \(i Q 2 (Bt + 2 Dt) - | V ft z(2) + 



^-^))^AtP , (61) 



where z(a) = X^^Li l/ n<1 is the Riemann Zeta function and the explicit expressions of the Gij have been given 
in sec. (2) of |3J. For a reference use in phenomenological analysis, let us also give the explicit expressions for the 
coefficients, which can be checked against the second and third order computations of the u spectrum as soon 
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as the latter will become available. We report only the coefficients differing from the starting ones, Gjjij ^ Gif 

C F z{2) 



G 



U22 



G\J2\ 



G 



U33 



G 



U32 



C F 


n f 


IT 2 


48 


C F 




IT 2 


n f 


Cf 


11 


7T 3 


~ 432 



2 

z(2) 
6 



C A [- 

vc A 

95 
216 







96 + 


f) 


^215 

vis" ~ 


17,2(2) 
12 


z(2) 
12 





4 

2471 
1728 



C F 

11^(2) 
24 



3 

32 



z{2) 



z(3) 



C f 



1 

16 



z{2) 



11 



C F C A z(2) + -C F z(3) 



40 



+ Cf nf 



_ £(2) 

' I 432 36 
/2971 17z(2 



z(3) 



23177 ( 155 z(2) llz(3) + 5z(4) 



C% 



10368 96 

7z(3) , z(4) 



V5184 48 12 
1541 4z(2) 77z(3) llz(4) 



12 



1728 



864 



9 



48 



16 



(62) 
(63) 



(64) 



1709 19z(2) z(3) 



72 



24 



(65) 



where CU 
obtained: 



3 and rif is the number of active flavors. Let us make a few comments about the results 



it is remarkable that we could explicitly compute Gu32 with the knowledge of the first two orders of the 
coefficient function Ch{w; a) only (cfr. eqs. (|22J) and (J23J ) - We can compute also the coefficients Gun,n-i 
for n > 3, i.e. the exponent Gjj to NNLO. As explained in detail in NNLO means indeed that for 

each power of a we can compute the three principal logarithms: 5 



Gun,n+1 Ct n L n , Gun,nCt n L n , Gu n ,n-1 Ot n L n . (67) 

By general counting arguments, one would expect that the NNLO corrections to also require the 
knowledge of the NNLO contribution to the coefficient function, i.e. of the 0{a 2 ) term Cg (w). That 

(2) 

is actually not the case because the NNLO contributions proportional to C H (w) cancel between the 
numerator and the denominator in the definition of the effective form factor in eq. (|50|) . However, a 
complete NNLO resummation of the u-spectrum also requires the knowledge of the second-order correction 



to the coefficient function C\j , and for that C H (w) is needed (see ea. (|49|l '): 

Since the coefficients of the threshold logarithms in the u distribution and in the photon spectrum in the 
radiative decay JSJ differ from two loops on in NLO, the cancellation of long-distance effects in the ratio 
considered in 0] 

r> ( \ dT R /dt s (t s = u) 

R{u) = m 7~j / \ 68) 

dl jduyu) 

is not exact but occurs only in leading order. As previously discussed, the u distribution has additional 
long distance effects with respect to the radiative decay related to small hadron energies. 



5 If we look instead at the form factor itself, Sy, NNLO means that we can compute for any n all the coefficients with 
- 1 < k < 2n: 



V 1 ~n j 2n ~ Tl r 2n-l „to jn \- „n r n — 1 

i(Jn,2n« L, , L,X7n,2n-l CX L , ■■■ , 2-iUn,n & Li , 2-*Un,n-l a Li 

For instance for n = 3 we can compute the coefficients of the logarithms from power six down to power two included. 



(66) 
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The differential spectrum in u is obtained from the event fraction Rjj(u) by differentiation: 

-— = Cu(a)au (u; a) + du(u; a), 
1 du 

where we have defined 



ou (u; a) = —Y,u{u;a) 
du 



djj{u\ a) 



du 



Djj(u; a). 



(69) 



(70) 



The coefficient function is clearly the same in the partially-integrated spectrum and in the differential one, while 
the remainder function is obtained by differentiation and reads: 



Cf 

7T 



36 - 8u - 8u 2 +u 3 7 - 56u - 29u 2 + 8m 3 



12 



70 



logtt 



(71) 



3 Hadron mass spectrum 



In this section we resum to NLO the distribution in the invariant hadron mass squared, i.e. the distribution in 
the variable 



t = 



'x 



(0<t<l), 



(72) 



< (1 + ^) 2 " 

where in the last member we have kept the leading term for u — > only. This distribution is obtained by 
integrating the distribution in the hadron variables u and w with the previous kinematic constraint: 



1 dT 
f~dt 



D 



dudw — ; — —5 
r dudw 



t 



u w 



(1 + uf 



the integration covers the whole phase space D of the hadron variables: 



dudw 



dw 



du 



du 



l+u 



dw. 



I max[0,w — 1] 

It is convenient to evaluate the event fraction, defined like in the previous distribution as: 



R T (t) 



df 



rdt> 



with the end-point values Rt{0) = and -Rt(I) = 1- Integrating both sides of eq. J73J), one obtains: 

„,.,,,2 i 



R T {t) 



, , 1 rf2r 

dudw — — — — ( 
D r dudw 



t - 



(73) 



(74) 



(75) 



(76) 



By inserting the resummed form for the double hadron distribution and neglecting at first the remainder 
function, we obtain: 



Rr(t; a) = I dudw Ch(w; a)a[u; a(wmi,)] 



D 



t 



(1 + u.y 



+ 0(t,a), 



(77) 



where by 0(t, a) we denote terms which are zero at the threshold and vanish for a — 0. In order to isolate the 
large logarithms, let us simplify the domain D into the unit square 



< u, w < 1 

and simplify the kinematic constraint as well as: 

6[t-uw 2 /(l + u) 2 ] -> 



9[t- 



(78) 
(79) 
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Let us observe that the variable t keeps unitary range even after such approximations. We then obtain: 



R T (t;a) 



o Jo 



dudw Ch(w; a) a [u; a{w mi,)] 9\t — uw 2 ] + 0(t, a) 



dwCH(w;a) + / dwC[i(w;a)'E[t/w 2 ; a(wrrib)] + 0(t,a). 



(80) 



Let us remark that the present case is more involved with respect to the one treated in the previous section, 
because the hadron mass squared t is a combination of both variables used for threshold resummation, namely 
u and w. Neglecting infinitesimal terms for t — > 0, the expression above can be further simplified by neglecting 
the first integral and integrating the second integrand over w down to zero, 



Rr(t; a) 



dwCH(w;a)T,[t/w 2 ; a(wmt)] + 0(a,t). 



(81) 



By inserting the first-order expressions for the coefficient function and the form factor given in the previous 
section, we obtain: 



dw Ch{w] a) £ [t/w 2 ; a(w rrib)] 



1 



a Cf 



1,2 31 , 637 

-log^i H logt H 

2 8 12 B 144 



+ 0(a 2 ). 



(82) 



As anticipated in the introduction, we have a different coefficient for the single logarithm with respect to the 
hadron mass distribution in the radiative decay (JSJ or the u distribution of the previous section. 



As with the u distribution, we introduce a resummed form of the event fraction as: 

R T (t; a) = C T {a)^ T (t; a) + D T (t; a). 
All the functions above have an expansion in powers of a: 



(83) 



C T (a) = 1 + aC ( T 1] + a 2 C (2) + 0(a 3 ); 
E T (u; a) = 1 + (t^'W +a 2 Z (2 \t) + 0(a 3 ); 
D T (t; a) = aD [ rp(t) + a 2 D (2) (t) + 0(a 3 ). 



(84) 
(85) 
(86) 



Let us consider a minimal factorization scheme, where only logarithms in t are factorized in the effective form 
factor Yjt- Since we will consider a different scheme later on in this section, let us denote the quantities in the 
minimal scheme with a bar. The first-order corrections to the form factor and the coefficient function in the 
minimal scheme read: 



££>(*) = 

7T 



1, 2 31 

- log t H log t 

2 6 12 S 



a 



(i) _ 



C F 637 
~V 144 



= - 1.87744. 



(87) 



Note that the correction to the coefficient function is very large: for a(rrib) = 0.22 it amounts to — 41.3%. 

We now expand the resummed result in powers of a and compare with the fixed order result, which is known 
to full order a [HE]: 



R T (t;a) = 1 + aR^\t) + a 2 R^ } (t) + 0(a 3 ), 



with 

R^\t) 



Cf 

IT 



- log 2 1 
2 B 



31 



lost + 
12 B 144 



637 97 . 



18 



t + 



25 
18 



61 
144 



+ 



t 4 lost 



(89) 
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We obtain for the remainder function in first order: 



-m, , 97 25 o 61 , /5 3 , 1 A , 

4 1} W ^ Is* " 18*^ + 144* " (l* "2^+6* T 0g " (90) 

Since 

E T (1; a) = 1, (91) 

taking t = 1 in eq. (|83|) . we obtain a relation between the coefficient function and the remainder function in the 
endpoint: 

C T (a) = 1 - D T (1; a). (92) 
It is a trivial matter to verify that the above relation holds true for our first-order expressions. 

A compact definition of the minimal scheme can be given as follows. Since 

dw Ch(w; a) E [t/w 2 ; a(w m&)] = Ct ip) St (t; a) , (93) 

the coefficient function in the minimal scheme can be defined taking t — 1 in the above equation and using (|91|l : 

C T (a) = [ dwC H {w;a)J:[l/w 2 ; a(wm b )] . (94) 
Jo 

The effective form factor then reads: 

^ dwCniw-^^^/w 2 - a{wm b )\ 

i^T [t; a) — — j . 

J du; Ch{w; a) S [1/w 2 ; a(w mj,)] 

Let us comment the above result. The effective form factor involves a convolution over the hadron energy 
w of the coefficient function and the universal form factor, which cannot be reduced to an ordinary product 
by the standard moment transform. That is because the variable w enters not only in the first argument of 
S = E \t/w 2 ; a(wmb)] but also in the scale of the coupling a = a(wm,b). Analogously to the u distribution, 
there are long-distance effects in the effective form factor related to small hadron energies w <C 1, which are 
suppressed by the coefficient function. In the present case, however, there is an additional mechanism suppressing 
the small energy contributions: since E is evaluated in t/w 2 , small w's correspond to a large argument u = t/w 2 
of E(u), where there are no large logarithms, log l/u ~ 0(1), and one is inclusive at the parton level. We may 
say that this spectrum is "double protected" from the non-perturbative long-distance effects related to small 
hadron energies. 

The systematic expansion of the form factor is easily obtained by writing as usual: 

E T = e° T , (96) 

one obtains: 

oo n+1 

G T (t; a) = Y, J2 ° Tnk a " L t ' ( 9? ) 

71=1 fe=l 

where 

L t = log j (98) 
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and 



Gt\2 
Gtw 

Gt23 

Gtii 
Gti\ 

Gt34 
Gt33 



g 



T32 



G\2\ 

Gn + |i4i; 
6 



— G23; 



G 22 + ^Aj + jUiA,; 



>>-£" 



23 7 
--^i/3o + ^1(^1 + ^1) + 



= G 2 i + -A 2 + A\ 

7T 

= G34; 

= G32 + 5 A 2 ft, + 5 A 2 D, + H V (B, + D, ) + .4! 2 ft 

216 



3 --^ u , 3 

+ v v ( w - 12 z(4) ) + ^ Al ^ 2 - i^ 1 ( 23 * a - 15 *) + ^ Al * ( Bl + 2Dl ) 



/35z(2) 47 
1^2 



27 



(99) 
(100) 
(101) 
(102) 



(103) 
(104) 
(105) 



, . C F /547 A ,„,\ , 3 / 1117 7z(2) 5z(3) 
+ ~ ( 108 ~ ( +Al ( _ 1296 12 fi 



(106) 



For comparison with future higher-order computations, let us give the explicit values of the coefficients Gnj 7^ 
Gn: 



Gtu 
Gti2 



G 



T21 



G 



T33 — 



G 



T32 — 



31 C F 

12 7T ' 
7T 2 



~^ + Cf \Ya-—) +Ca l32 + — 



7T 2 

cv 

7T 3 



z ( 2 ) 83 \ ^ /941 11 7 , o ,\ ^ /107 11 «(3) 



144 



ni + C f Ca\ 

1296 f V 144 



288 



ll^(2)V r „ ( 5 , «(2) 



32 6 
+ 



(107) 
(108) 

5(109) 



,25 «(2) 
^^^-162"^ 



+ Gi 



1057 llz(2) \ 2 /_83_ £(3) 
5184 24 7 F V 648 3 



(110) 



Cf 

7T 3 



263 z(2) 



2592 36 



n f + CFCA {wm + 



( 28493 457^(2) 77^(3) 5 2(4) 



96 



12 



+ 



+ 



f<i n 2353 _ 47z(2) + llz(3) ^ + ^ f 60287 + 7z(2) _ 31z(3) _ 47z(4) ^ 



5184 



48 



12 



V 5184 12 



12 



+ C A n s [- 



6515 17z(2) z(3) 



5184 



36 



24 



/ 31841 229z(2) 77z(3) 11 z (4) 



V 10368 



144 



48 



16 



.(111) 



Let us make a few remarks: 



• The coefficient of the single logarithm at 0(a) is different from the previous case as well as from the 
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radiative decay: 

Gtu 7^ Gun — Gui (112) 

because Gtu takes a kinematic contribution from A±, i.e. from the double logarithm at one loop in E. 

The logarithmic structure of Gt radically differs from that of G because the integration variable w enters 
not only the argument of the running coupling a = a(w rrib) but also the argument of the logarithm 
L = log(t/w 2 ); 

• If the hard scale was set by the heavy flavor mass, Q = nib ( a kind of frozen coupling case with respect to 
the real case) we would have the following values for the coefficients: G^ 22 = G22 + 7/24 A\ + 5/4 A\Pq 
and G T 33 = G 33 + 83/ 648 A\ + 7/8 A\[3 a + 35/18 A^ 2 . 



As discussed in the previous section, in phenomenological studies one may want to replace the perturbative 
expression of S(u; a) with a fit to some experimental data or with a non-perturbative model. The minimal 
scheme cannot be used directly in these circumstances because the effective form factor a) involves the 
integration of S(u; a) in the unphysical region u > 1. There are various ways to deal with this problem. One 
way could be for example replacing S(it; a) with the non-perturbative quantity T, np (u; w) for u < 1, while still 
keeping the perturbative S(it; a) in the unphysical region u > 1. The perturbative form factor is indeed an 
analytic function of u, which can be continued to any value of u. To avoid the inclusion of the perturbative 
S(u; a) for u > 1, let us consider instead a non-minimal scheme with an effective form factor defined as: 



S T (*; a) 



fydwCniw; a)S[i/w 2 ; a(wmb)] 
dw Ch(w, a) 



(113) 



where the standard form factor S(n; a) has been extended to arguments larger than one, u > 1, since: 



a) 



S(m; a) for u < 1; 
1 for u > 1. 



(114) 



Because of the definition, it holds: 



Explicitly one has (cfr. the r.h.s. of eq. I|8()|l): 



E T (1; a) = 1. 



(115) 



1 pVi A 

dw Ch(w; a)!] [t/w 2 ; a(w nib)] = / dwCff(w;a) + / dw Ch(w; a)T, [t/w 2 ; a(w nib)] ■ (H6) 



The coefficient function is given in the new scheme by: 



C T (a) 



dw Ch(w, a). 



(117) 



By inserting in eq. (|113f) the perturbative expansions for Cr and for S, we obtain: 



4 1} (*) 



1,2 31 

- log t log t - 

2 B 12 B 



151 - 232 1 3 / 2 + 81 i 2 



72 



(118) 



and 



r iX) _ — 



C F 335 



7T 144 



0.98735. 



(119) 



The first of the above equations shows that St, unlike St, is not defined in a minimal factorization scheme, 
i.e. it does not contain only logarithmic terms a n L\, but also contributions of a different form. Note that 
the coefficient function in (|119fl has a much smaller value than in the minimal scheme, giving a hint of better 
convergence of the perturbative series in the modified scheme. 
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Matching with the first-order spectrum, one obtains for the remainder function in the modified scheme: 

(120) 



7T 



97 29 0,0 9 2 25 o 61 , f5 3 , 1 A \ , 

— * i 3/2 + -< 2 < 3 H t 4 - -f - -t 2 + -t 4 logi 

18 9 8 18 144 \3 2 6 f 



The coefficients of the logarithms in the exponent of the form factor Gt are the same in the minimal scheme 
and in the modified one: 

G Ti j = G Tij for j > 1. (121) 

The non logarithmic coefficients in the modified scheme are given by (Gtw = and Gt20 — by definition of 
minimal scheme): 

23 5 

Gtw = A x + -(B 1 + D 1 ); (122) 

36 6 

G T2 o = -§A a + ^r(Bi+ D 1 ) 2 + Ul3 2 + D 2 ) + ^:f3 B 1 + 
36 24 6 36 

47 ( 5z(2)^ ( ^ /2057 23z(2) 5z(3) 



+ A 1 (B 1 + D 1 )[ + — —l+Ai, 

K J \ 54 6 J 1 V2592 36 6 

Cf ,„ ^ s /547 A ,„,\ C F . / 90121 102(3) 
+ v (5 1+ ^)( M -4, ( 3)j + -A 1 ^^ r + -li + 12,(4)). (123) 

Explicitly: 

G T10 = (124) 

- 5{ c ,(-|! + 5 =(2) + |,3,) + ,,(i|i-|^ 
♦*(-^-> + £«»> + - <*>)}■ 

The relations between the coefficient function in the minimal scheme and in the modified one are obtained 
imposing that 

C T e GT = C T e° T + 0(t;a) (126) 

and read: 

fi\X) ft 



T — °T — ^TIO: 

ft{^) fi{2) ftW ft ~ ft'A ft /'107^ 



2 

The first equation can be directly verified by inserting our first-order expressions. The second order correction 
to the coefficient function is unknown at present in either scheme and its determination requires a full two-loop 
calculation; the second of the above equations simply allows us to transform the coefficient function from one 
scheme to another one. 

The Babar collaboration has recently presented the differential spectrum in t which is obtained from the 
previous one by differentiation: 

1 dT d „ , . 

vTt = Jt RT ^ (128) 

Due to their relevance, let us present explicit formulas. The resummed spectrum reads: 

i c/r 

F -77 = C T {a)a T (t; a) + d T (t; a) (129) 
1 dt 

where 
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• the coefficient function is the same as in the event fraction, since it is independent on t: 

• the effective form factor is 



ar(t; a) = —Y, T (t; a) 



More explicitly: 



a T (t; a) = 



j\j:dw/w 2 Ch(w; a)a[t/w 2 ; a(wmb)\ 
J dwC H {w\ a) 
du/(2\/Tu) Ch y\/t/u] a(wrrib)j a 



Jg 1 dw Ch(w; a) 



where a double representation as an integral over w or over u respectively has been given; 
the remainder function is: 

d T {t- a) = -rD T (t; a) = ad^ ] (t) + a 2 d^ (t) + 0(a 3 ). 

The explicit value of the first-order correction is: 



4 1} (<) = ^ 

7T 



67 29 r 15 25 , 55 , 5 , , 2 

\ft H 1 1 2 H r log t + 3t log t 1 3 log t 

18 6 4 6 36 3 B 5 3 6 



(130) 



(131) 



(132) 



(133) 



Let us make a comment about the non-perturbative effects entering the differential mass distribution. The 
expression for the effective form factor ot(^; ol) at the second member of eq. <|131fl involves an integration over 
w from y/i up to one. Since the running coupling is evaluated in Q = w nib, the smallest hard scale contributing 
to the distribution is 

Qmin = Vi nib = m X - (134) 

In order to avoid the infrared pole in the coupling — the well-known Landau pole — implying a breakdown of 
the perturbative scheme, one has to impose the condition 



(135) 



which is also very reasonable from the physical viewpoint. Resummed perturbation theory therefore signals that 
the hadron mass distribution cannot be computed for hadron masses of the order of the hadron scale because 
of the appearance of the Landau pole 6 . 



4 Conclusions 

In this work we have presented next-to-leading resummed expressions for the distribution in the final hadron 
mass/energy ratio and for the distribution in the invariant hadron mass in the semileptonic decays 

B -> X u + I + v. (136) 

By expanding our formulas, we have obtained the coefficients of all the infrared logarithms to 0{a 2 ) and of the 
leading ones to 0(a 3 ). These two spectra have different logarithmic structures from each other, which are both 

6 The universal QCD form factor a [u; a(Q)] has an infrared singularity at u = exp [— 1/(2 /3q a(Q))] , related to the Landau pole 
for mx ~ \/ AqcdQ S> ^qcd- We assume that the latter has been regulated in some way; for a recent discussion see for example 

nn . 
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different also from that one in the radiative decay (JSJ. That occurs because these spectra involve integration 
over the total hadron energy Ex, which sets the hard scale Q of the hadronic subprocess in (|136[) : 

Q = 2E X . (137) 

Long distance effects manifest in perturbation theory in the form of large infrared logarithms, which are usually 
factorized and resummed in QCD form factors. Universality of long-distance effects therefore shows up in 
perturbation theory as the occurrence of equal form factors in different distributions. That implies that the 
spectra we have considered in this work have different long-distance effects from each other as well as from the 
radiative decay (JSJ. There is no simple connection between these semileptonic spectra and the hadron mass 
distribution in radiative decays. For both spectra, we have introduced effective, i.e. process dependent, form 
factors, which factorize the large logarithms to all orders in perturbation theory. These effective form factors 
can also be computed in a phenomenological way by inserting, in place of the perturbative QCD form factor 
E [it; a(wrrib)], the form factor E np (it; w) computed with a non-perturbative model or a fit to experimental 
data. In the case of the distribution in the hadron mass squared t, we have also considered a non-minimal 
factorization-resummation scheme, which seems to have better convergence properties of the perturbative series 
with respect to the minimal one. 

There are other important semileptonic spectra which have similar properties to those of the distributions 
considered here, i.e. long distance effects which cannot be factorized into a process independent form factor, to 
be extracted for example from the radiative decay 

Finally, we have also shown that the cancellation of long-distance effects in the ratio constructed in pi] occurs 
only in leading order while it is violated at the level of a 2 log 2 (i) terms. 
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